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Abstract 

An apparent violation of hyperscaling at the endpoint of the critical line in the 
Domany-Kinzel stochastic cellular automaton finds an elementary resolution upon not- 
ing that the order parameter is discontinuous at this point. We derive a hyperscaling 
relation for such transitions and discuss applications to related examples. 
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Domany and Kinzel's stochastic cellular automaton (DKCA) continues to attract interest 
as a particle system affording a test of universality in nonequilibrium critical phenomena. [I], 0] 
In the phase diagram of the DKCA, a line of critical points separates the active phase from the 
absorbing vacuum phase. There is good numerical evidence |^, M that the critical behavior 
along this line is that of directed percolation (DP), except at one terminal point (called com- 
pact directed percolation or CDP), where the asymptotic behavior is known exactly. Essam 
showed that the critical exponents governing the moments of the cluster-size distribution in 
CDP violate the usual hyperscaling relation, and derived the appropriate generalization. H 
In this note we point out that the hyperscaling relation amongst the exponents describing 
spreading also requires modification. This leads us to formulate a hyperscaling relation for 
discontinuous transitions into an absorbing state. 

In the one-dimensional DKCA, the state o~i(t + l) (=0 or 1), of site i at time t + 1 depends 
upon 7] = (Ti(t) + 0i_i(t). If we let Pr[oi(t + 1) = 1] = h(Ti), then h(0) = 0, h(l) = p u and 
h(2) = p 2 . (All sites are updated simultaneously in this discrete-time process.) For each 
p 2 < 1 there is a critical value pi )C , such that 6^ (the vacuum — all sites in state 0), is the 
unique steady state for p\ < pi jC , whilst for p\ > px sC there is also an "active" stationary 
state with a nonzero density of sites in state 1. It is believed that the critical behavior along 
the line Pi )C (p2) is the same as that of DP, for p 2 < 1. 

The case p 2 = 1, (CDP) is special. On this line, 5±, the all-1 state, is absorbing, as 
well as 5$. The evolution is determined by the coalescence of random walks marking the 
boundaries between strings of 0's and l's; for p\ = 1/2 the walks are unbiased. In fact, 
(pi = 1/2, p 2 = 1) is a critical endpoint: for p 1 < 1/2, 5^ is the globally attractive stationary 
state, and conversely for S\, when pi > 1/2. (By "globally attractive" we mean that for 
Pi < 1/2, initial configurations with an infinite number of 0's evolve to 5$ with probability 
1, and that there is a nonzero probability to reach 5^ if there is at least one site in state 
initially.) Let p\ = 1/2 + A. It is easy to see that CDP is invariant under the transformation 
which takes A to —A and interchanges 0's and l's. 
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Consider the evolution at (1/2,1), starting with the origin in state 1, all other sites zero. 
The number n(t) of l's at time t is a simple random walk on Z starting at n(0) = 1, with 
an absorbing boundary at the origin. From well-known properties of random walks, we have 
that the survival probability P(t) = Pr[n(t r ) > 0, t' < t) ~ t -1 / 2 , and that the mean-square 
population of l's in surviving trials is n surv (t) = E[n(t) 2 \n(t') > 0, t' < t] ~ t. It follows 
that the mean number of l's (averaged over all trials, including those which die out) is 
asymptotically 0(1). Spreading of a critical process from a localized source is customarily 
described by a set of exponents defined via the relations || 

P(t) oc t~ 5 , (1) 

n(t) oc f, (2) 

and 

R\t) oc t\ (3) 

where n(t) denotes the mean population size and R 2 (t) the mean-square spread of particles 
(l's, in our notation) about the origin. Along the critical line in the DKCA, the exponents S, 
r), and z take universal values (i.e., those of 1+1-dimensional DP), except at (1/2, 1) where 
evidently S — 1/2, T] = 0, and z — 1. 

The spreading exponents are connected by a hyperscaling relation || 

45 + 2r] = dz, (4) 

where d is the number of spatial dimensions. While the exponents for DP are consistent with 
this relation, those for the point (1/2,1) (with d = 1) aren't. Does this signal a violation of 
hyperscaling in CDP? Here it is useful to note that Eq (U) is a special case of a more general 
relation 
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where (3 is the usual order-parameter exponent, defined, for the DKCA, through pi oc (pi — 
Pi, c ) , for Pi > Pi,c] Pi is the stationary density of l's. f3' governs the ultimate survival 
probability (starting from a localized source): P^ = lim^oo P(t) oc (pi — Pi^ ■ Eq (||) 
reduces to Eq(^) when the order parameter and P^ are governed by the same exponent, as 
they are in the contact process and other DP-like models. Models with multiple absorbing 
configurations 0, and branching annihilating random walks with even parity ||, can have 
(3 7^ (3', in which case the exponents violate Eq (ffl), but satisfy Eq 



It is known that j3' = 1 in CDP. |], §, H |10| The order-parameter exponent, (3, by 
contrast, is zero. This is because (1/2,1) marks a discontinuous transition, as is readily seen 
by recalling the symmetry property noted above. Since p\ = for p± < 1/2, (the globally 
attractive state is 5^), it follows that pi = 1 for pi > 1/2. Strictly speaking, f3 is not defined 
here. But since lim/3\o A 13 = 1 for any A > 0, it is natural to associate the value (3 = with 
the discontinuous transition. Indeed, the values 5 = 1/2, z = d = 1, and rj = (3 = yield an 
identity when inserted in Eq (|5|). 

In fact, we can eliminate any reference to the ill-defined exponent (3 by adapting Grass- 
berger and de la Torre's scaling argument to discontinuous transitions. Consider a model 
with a transition from an absorbing to an active state at A = 0, with exponents 5, 7], z, and 
f3' defined as above. Suppose, however, that the order parameter p is discontinuous, being 
zero for A < 0, and 

p = p + /(A), (6) 

for A > 0, where po > 0, and / is continuous and vanishes at A = 0. The scaling hypothesis 
for spreading from a source postulates the existence of two scaling functions, defined via M 

p(x, t) ~ t v ~ dz/2 G(x 2 /t z , At 1 /"!! ), (7) 

and 

P{t) ~r s $(At 1/v "). (8) 
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(Here p(x, t) is the local order-parameter density. v\\ governs the divergence of the correlation 
time r at the transition: r ~ A~ v w.) Existence of the limit implies that ~ x@ as 
x — > oo, with f}' = 5v\\. In a surviving trial, the local density must approach the stationary 
density p as t — > oo, so p(a;, £) ~ A^'po, for t — ► oo with fixed x, and A small but positive. It 
follows that G(0, y) ~ y^' for large y. On the other hand, we must have G(0, y) ~ y- u \\( r i~ dz / 2 ) 
for lim^oo p(x, t) to exist. Comparing these asymptotic behaviors, we find a hyperscaling 
relation for transitions at which the order parameter is discontinuous 



5 + r7 = y . (9) 

The interpretation is immediate: simply note that 5 + rj is the exponent governing the mean 
population in surviving trials, and that the radius R t of such a cluster grows ~ t z l 2 . (Since 
the density is positive, clusters are compact, not fractal as in DP.) Eq (||) is nothing more 
than the scaling law for the volume of a (^-dimensional sphere of radius Rt . It should apply 
whenever power-law growth produces compact colonies. (By "colony" we mean the set of 
particles, or l's — in general, active sites — at time t, descended from a single active site 
at t — 0. "Compact" means that the density in surviving colonies remains finite as t —>■ oo.) 



A closely related example is the voter model. || |TT] In this continuous-time Markov process 
sites of Z d are either in state or state 1. If site i is in state 0, it changes to 1 at rate Ari, 
where r± is the number of nearest neighbors in state 1. (Note that only sites at cluster 
boundaries can change state.) Similarly, sites in state 1 change to at rate r . (The case 
with A 7^ 1 is usually called the biased voter model, for obvious reasons.) Clearly, 5^ (5\) is 
the attractive stationary state for A < 1, (A > 1), so p% jumps from to 1 at A = 1. The one- 
dimensional voter model is essentially a continuous-time version of CDP with p\ = A/(l + A). 
Thus we expect all results for critical exponents in the DKCA to apply as well to the voter 
model in one dimension. 

We can analyze spreading in the (unbiased) voter model in d > 2 as follows. Let n s 
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be the number of l's after s changes in the configuration; no = 1. n s is again a simple 
random walk on Z, with an absorbing boundary at the origin, so we know that the survival 
probability ~ s -1 / 2 , the mean-square displacement over surviving walks ~ s, and the mean 
displacement over all walks is 0(1). The latter implies rj = 0. The number of steps per unit 
time is proportional to the boundary (number of 0-1 nearest neighbor pairs), b s \ 



In one dimension, b s = 2, so t ~ s. For d > 2 the boundary depends on the shape and 
internal structure of a colony; we assume b s ~ n] as s — > oo. Then since n s ~ s 1 / 2 , we have 
t(s) ~ s^ 2 , n SMTO (t) ~ t 1 /(2-i) ) anc i p( t ) _ r V(2-7) > so that 5 = 1/(2 - 7). Since b < qn 
on a lattice with coordination number q, 7 < 1, implying that | < 5 < 1. 

We expect asymptotic properties to be captured by a continuum description, which for 
the voter model takes the form of a very simple stochastic partial differential equation, 

dp(x, t) 



dt 



= V 2 p(x,t)+r ] (x,t), (11) 



where p(x,t) is a coarse-grained density of l's and the (Gaussian) noise has autocorrelation 
< r](x, t)r](x', t') >= 5(x — x')5(t — t')p(x, t)[l — p(x, t)]. (There are no fluctuations when the 
density is pinned at one of the absorbing values.) In mean-field approximation (neglecting 
the noise) the population spreads diffusively and z — 1. But we know that z = 1 even for 
d — 1, suggesting the mean-field exponent is correct for all d. 

Simulations of the voter model in two and three dimensions support this conjecture. In 2- 
d, a study of 2 x 10 6 independent realizations (all starting with a single 1), up to a maximum 
time of 10 4 , yield 5 = 0.95(3), rj = —0.01(1), and z = 0.99(2) (figures in parentheses indicate 
statistical uncertainties). In three dimensions (10 6 realizations up to t = 900), we obtain 
5 ~ 0.96(2), i] = 0.00(1), and z = 1.00(1). (In both two and three dimensions the colony 
mean-square radius of gyration grows linearly with time.) Moreover, the boundary exponent 
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7 appears to be unity in each case. In 2-d, 7 shows considerable time dependence, increasing 
monotonically from 0.82 for t = 1000 to 0.93 for t = 20, 000, the longest study attempted. 
In three dimensions our study of colonies surviving to t = 900 yields 7 = 0.975(4). The 
simplest conclusion is that 7 = 5 = 1 for d > 2. Then the hyperscaling relation, Eq([|) is 
violated for d > 2, marking d = 2 as a kind of upper critical dimension for spreading in voter 
model. (Consistent as well with the slow approach to asymptotic scaling found in the 2-d 
simulations.) The density within a colony ~ n surv (t) / Rf ~ f}~ d l 2 when r\ = and 5 = z = 1, 
so for d > 2 colonies are not compact, and Eq(|9|) is not expected to hold. 

In summary, we have shown how hyperscaling applies to critical spreading at discon- 
tinuous transitions, and discussed applications to compact directed percolation, and the 
closely-related voter model. 

We are grateful to N. Inui, M. Katori, and H. Takayasu for helpful discussions. 
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